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Abstract — To achieve high range resolution profile (HRRP), the geometric theory of diffraction (GTD) 
parametric model is widely used in stepped-frequency radar system. In the paper, a fast synthetic range 
profile algorithm called orthogonal matching pursuit with sensing dictionary (OMP-SD), is proposed. It 
formulates the traditional HRRP synthetic to be a sparse approximation problem over redundant dictionary. 
As it employs a priori information that targets are sparsely distributed in the range space, the synthetic 
range profile (SRP) can be accomplished even in presence of data lost. Besides, the computational 
__i complexity is reduced by introducing sensing dictionary (SD) and it mitigates the model mismatch at the 

^vq same time. The computation complexity decreases from O(MNDK) flops for OMP to 0(M(N + D)K) 

^ flops for OMP-SD. Simulation experiments illustrate its advantages both in additive white Gaussian noise 

(AWGN) and noiseless situation, respectively. 
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1. Introduction 

A HRRP is the phasor sum of the time returns from different scatterers on the target located within 
a resolution cell. From a geometric point of view, a HRRP represents the projection (in range) of the 
apparent scattering centers onto the radar line of sight (LoS). It has been used into radar signal processing 
__-__. In GDT model, the radar target is no longer a point but composed of multiple scatterers along with 
radar LoS. Different targets are characterized by different scattering mechanisms and they are expressed 
as a function of frequency. From the perspective of digital signal processing, a process to identify radar 
target in GTD model is just the same as the process to estimate model parameters (containing scattering 
mechanisms, intensity and scatter range cells). In other words, it is also a process to estimate parameters 
of position and magnitude of target scatterers in range gate. While, in realistic process of synthesizing 



January 27, 2013 



DRAFT 



DRAFT 



2 



range profile, the returned signals are always inevitablely interfered by passive or/and active jamming. 
For the case, many returned signals are corrupted or even invalid. If signal sample is invalid, it has to 
be discarded. If exploiting the measured echo signals to synthesize range profile with traditional SRP 
algorithms directly, such as O-Q, the SRP of target is either incorrect or flashed from coherence 
processing internal (CPI) to another one. Although interpolation or extrapolation strategy are helpful 
to improve the SRP, the range resolution remains constrained by system band width. However, it's 
worth noting that the significant physical scatterers are sparse in actual targets, which implies that strong 
scattering cells are also sparse for the target's SRP JSJ. This is consistent with sparse signal representation 
of compressed sensing (CS) theory appeared in recent years [9]-|13|. Hence, in the signal representation 
point of view, synthesizing range profile is equivalent to recovering a high-dimensional sparse signal form 
a low-dimensional measurements, usually accompanied by samples loss. It is an undetermined system. 
This will be confirmed in next sections. 

In the past few years, some algorithms have been developed to solve an undetermined system by using 
sparse property. They are generally grouped into two categories in CS community, i) minimum l\ norm 
reconstruction, i.e., optimization based on the l\ norm can exactly recover sparse signals and closely 
approximate HRRP with high probability. This is a convex problem that conveniently reduces to a linear 
program known as Basis Pursuit (BP) 0; ii) greedy algorithms, one representation of those is orthogonal 
match pursuit (OMP) lfT4l - |[T5ll : Considering it has a substantial gap between the computational cost of 
OMP and the cost of BP, we develop faster reconstruct algorithms than OMP algorithm in the paper. For 
describing convenience, the measurement matrix in measurement system (see Eq. (3)) is called dictionary. 
Each column in dictionary is called an atom. Meanwhile, it calls that it is if-sparse if containing K 
nonzero entries in a vector. 

However, in the GTD model with multiple scatterers, a few scattering mechanisms should be considered. 
With the increase of atom number in dictionary, the computational cost increases. Although a simplified 
scattering model can be used to approximate multiple scatterers model as discussed in section 3, model 
mismatch can degrade the success recover probability, which deteriorates the cumulative distribute error 
(CDE) of SRP. Similar to the SD in |[T6l - |[T7l which is used to mitigate inter-atom interference (IAI), the 
SD is introduced to mitigate model mismatch in this paper. To the authors' knowledge, in the previous 
work, there is not report in GDT model yet. Using SD, it can reduce computational complexity and 
mitigate the model mismatch so as to improve the recover probability of SRP. The main contributions 
of the paper are three aspects. Firstly, for measurement data loss, it adapts sparse property of HRRP 
to synthesize range profile. In the second, it mitigates model mismatch by introducing SD. Thirdly, an 
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improved fast algorithm (i.e., OMP-SD) is proposed. 

The paper is organized as follows. In section 2, it first presents the GTD scattered model in frequency 
domain and then, establishes measurement system in stepped frequency radar (SFR). After that, it briefly 
reviews existing algorithms to solve the model and presents approximate OMP algorithm (A-OMP) in 
section 3. In section 4, it presents a strategy to construct SD. It mitigates model mismatch effectively. 
Besides, a fast algorithm (OMP-SD) to synthesize HRRP is proposed. Monte Carlo simulations illustrate 
the performance of the proposed algorithm both in AWGN and noiseless situation respectively, in section 
5. Finally, some conclusions and further work are provided in Section 6. 

Notation: It denotes vectors and matrices by boldface lowercase and uppercase letters, respectively. 
Uppercase Greek letters also represent matrix in this paper. (-) r denotes the transpose operation, (•)•" 
denotes the conjugate transpose operation, Further, || • \\% refers to the 1% norm for vectors. || • ||oo refers 
to the Zoo norm for vectors. The vec(-) operator vectorizes a matrix by stacking its columns. R € R LxA/ 
and R G C LxM denote a real- valued and complex- valued matrix and let 5t{-} and be real part and 
imaginary, respectively. (•)+ denotes the M-P generalized inverse. 

2. Problem Formulation 

In this section, it briefly presents the GTD scatter model of SFR return signal. SF pulse trains are 
created by transmitting a train of M identical baseband pulse with different carrier frequencies. The 
carrier frequency of the m-th (m = 0, 1, • • • , M — 1) pulse is f m = fo + A/, where /o is the initial 
frequency and Af is the frequency step size. In the stretch processing EOl . ll2~Tt . the range resolution 
is Ar = c/(2MA/), and the ambiguous range AR = c/(2A/) (c is the speed of light). The two- 
dimensional geometry of the radar scenario is shown in Fig. 1. For the convenience of signal modeling 
and derivation, it is assumed that the target is stationary and it falls in the range gate [L,L + Lq] in one 
CPI, where, L = QAR, and Lq = NAr). (Q and N are nonnegative integers). Meanwhile, it assumes 
that the target can be present only the grid points and let us discretize the range space by Ar in Lq. 



In one aspect angle, the parametric GTD scatter model of SFR at frequency f m can be represented as 

follows m-iEa, ea, 




d=l 1=1 



(1) 



where, 




(2b) 



(2a) 
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Fig. 1: Scenario of a target in SFR 



In (1), Gd, oid are the complex amplitudes, scattering mechanism of <i-th scatterer, respectively. rai 
denotes range space w.r.t. /-th range resolution cell in m-th pulse. D is the number of scatterers. u m is 
the AWGN with mean zero and variance a 2 . 

The y m in (1) consists of N uniformly sampled time-domain data from the baseband echo signal of 
m-th pulse (N = 2Lo/(cAt) and At = 1/(MA/)). The model can be written into a matrix form as 
follows, 

y = $x + u (3) 

where, y G C Mxl , $ G £MxDN and x g ^DNxi are measurement vector, dictionary and HRRP 
index of the target, respectively. <J> = $2, • • • , $£>] and 

[*d]m,n = G d [j(l + mA/// )p • exp(-j2^/ m (r + n/ (MA/))) (4) 

x = [x^f x^ ••• x^] T , Xrf G M Afxl , arf £ (fi is a set composed of scattering mechanisms), 
u G C Mxl is the AWGN vector. For convenience of the later describing, it defines $ = [<3?i, $2, • ■ , 3>d] 
and $>d — [<j>dli 4'd2, ■ • • i^dw]- denotes <i-th block matrix of denotes the i-th atom of &d- 

d G A = {1,2, • • • ,-D}. All atoms are normalized throughout the paper. r$ is the radial distance from 
radar antenna to reference point on the target. In realistic settings, M « N < DN, hence it is an 
underdetermined system in (3). It is to reconstruct of a high-dimension sparse vector x from a small 
number of linear measurements y and dictionary <£. 

3. The A-OMP Algorithm 

To solve an underdetermined system of linear equations in the above form (3), recently Zi-norm 
minimization as an effective technique has attracted attention in the CS community ||9l- |fT3l . It is a 



January 27, 2013 



DRAFT 



DRAFT 



5 



convex optimization problem in noise setting: 

(Pi :) min ||x||i subject to ||y — <3?x||2 < e (5) 

X 

For the problem Pi, it can be solved by linear programming (LP). Nevertheless, general-purpose LP 
solvers require about 0(D 3 N 3 ) flops. Thus, it is beyond the scope of discussion in the paper. As a matter 
of fact, many of the applications of Pi can be attacked heuristically by fitting sparse models, using greedy 
stepwise least squares. A widely used algorithm for sparse signal recovery is the OMP algorithm for the 
recovery of the support of the K-sparse signal in (3), which requires O(DNMK) flops llT4l - |[l"5l . For 
noiseless case, the exactly recovery condition (ERC) of OMP was derived by Troop in 2004. T.cai et al. 
derived a new ERC both in the bounded noise and Gaussian noise ||23l . 

In (3), considering the computational cost, the sparse solution can be solved with OMP algorithm 
directly. Rather than minimizing an objective function in (5), OMP constructs a sparse solution to a 
given problem by iteratively building up an approximation, the vector y is approximated as a linear 
combination of a few atoms in dictionary <£, where the activeset of atoms to be used is built column by 
column, in a greedy fashion. At each iteration, a new atom that best correlates with the current residual 
is added to the activeset. The standard OMP algorithm can be found in lfT4l . 

For GTD multi-scattering model, a few scatter mechanisms are considered. It increases the atom number 
in dictionary, and hence it increases computation. To decrease computation caused by multiple scatterers, 
it's a straight way to synthesize range profile of target that using a single scattering mechanism instead 
of multiple scattering mechanisms (i.e., to replace with <£d)- 

Just as the description in section 2, the atoms in dictionary (i.e., columns of $<j) are normalized so 
that ||<^di||2 = 1, for i = 1, 2, • • • , N. It denotes by c C S = {1, 2, • • • , N} the support of x^, which 
is defined as the set of indices corresponding to the nonzero components of x^. For matrix <3?d, &d(c) 
denotes the matrix formed by picking the atoms of <I>^ corresponding to indices in set c. Following 
the same convention as section 2, 4>di represents the i-th atom of It calls (pdi a correct atom if the 
corresponding x^ ^ and call 4>di an incorrect atom otherwise. With slight abuse of notation, we use 
$rf(c) to denote both the subset of atoms of <3?d with indices and the corresponding block matrix of <J>^. 
A detailed description of approximation orthogonal matching pursuit (A-OMP) algorithm is presented as 
follows. 

Similar to OMP, the A-OMP is a stepwise forward selection algorithm and is easy to implement. A 
key component of A-OMP is the stopping rule which depends on the noise structure. In the noiseless 
case the natural stopping rule is r« = 0. That is, the algorithm stops whenever r« = is achieved. In this 
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Algorithm 1 :A-OMP 
Input: 

The measurement vector, y; 
The dictionary, d € A; 
the error threshold, e; 
Main Procedures: 

1: Initialize the residual ro = y and initialize the subscript set of selected atom Co is empty. Set i = 1. 
2: Find the atom 4> t . that solves the maximization problem 

ti = max |</>^rj_i|, (t G 5, 0^ is t/ie ii/i atom m 

and update q = Cj_i U {^}. 
3: Let Pj = $rf(cj)(<I> ( i(Q) i ^<I> ( i(cj)) _1( I>rf(cj)^. Denote the projection onto the linear space spanned by 

the elements of $^(ci). Update r.j = (I — Pj)y. 
4: If the stopping condition is achieved (i.e., ||r-j||2 < e), go to 5. Otherwise, set i = i + 1 and return 

to 2. 

5: Pick out the range scattering cells w.r.t set q. 

6: Calculate the scattering intensity in these range cells determined in the previous step with P^y. 
7: Reconstruct SRP using the scattering intensity and range scattering cells. 
8: Return SRP. 



paper, both noiseless and the case of AWGN in which Uj ~ A(0,cr 2 ) are considered. The stopping rule 
for each case and the properties of the resulting procedure are discussed in article |[T4ll . 

As a special case of multi- scattering center, for a single scattering mechanism, the A-OMP algorithm 
procedure is the same as OMP, but it has a significant different physical meaning. Because the dictionary 
<l>d just as a sub-block of $ in (3). Thus, it is called approximate OMP in the paper. Once the subscribe 
set is determined with A-OMP or OMP, the SRP can be obtained with Least-Square (LS) solution. 

4. OMP Algorithm via Sensing Dictionary (OMP-SD) 

In (3), the sparse vector x can be obtained by OMP or A-OMP directly. However, two major problems 
cannot be avoided in this case. For OMP, it has to search all atoms in dictionary $ to find the best matched 
atom at each iteration (in the paper, the dictionary is M-by-DN dimension matrix); For A-OMP, it just 
needs to find the best matched atom in M-by-A^ dimension dictionary at each iteration, but it leads to 
model mismatch and increases CDE of SRP. Thus, an improved algorithm via sensing dictionary (i.e., 
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OMP-SD) is developed to overcome drawbacks of both OMP and A-OMP There are two advantages 
with OMP-SD to synthesize HRRP On the one hand, as a result of the atoms in SD are independent on 
scattering mechanisms, it mitigates model mismatch. On the other hand, it reduces computation because 
of the searching dimensional of dictionary reduced from M-by-DN down to M-by-AT. 
4.1. Dictionary Pre-processing 

For the convenience of following analysis, the dictionary of in (3) are divided into D M-by-N 
dimensional block matrix firstly, which are denoted by $i, $2> • • ■ > $D an d eacn block matrix d G A, 
corresponds to a different scattering mechanism. The Eq. (3) can also be rewritten as, 

y = [$i|$ 2 |---| $ o]x + u (6) 

In (6), it considers $i (w.r.t the 1st scattering mechanism) as an example. In ideal condition, the Gram 
matrix <&i<&d = I> but it is not the case because the dictionary is over complete, so it has to make 
$>i$>d ->• I, (d € A and d ^ 1), extremely, which needs to solve the problem max^g^ ||I — $f* 3><i||oo- 
According to the idea, it should find an M-by-N SD W ( being the same dimensional as block matrix 
<J?d ), which is independent on scattering mechanisms. The SD can be found by solving the problem P2 
of the follows, 



( min b\ + 7&2 (7a) 
w 

s.t p-dio^CW^Hoo <h (7b) 

s-t. ||p||oo < b 2 (7c) 



(ft 



vec((W H $ d ) k j),k£l (7d) 



KdeA (7e) 
In (7a), 7 is the regular factor. It sets 0.5 in the paper. Both 61 and 62 are unknown but determined 

variables. They reflect the IAI level between W and dictionary <J>^. As the problem P2 is a convex 

problem, the sensing dictionary W can be obtained offline with efficient algorithms. There are many 

software pockets to solve the problem such as cvx [24 1 etc. 

In P2, for the first constraint, it means maximizing correlation of diagonal elements of matrix which is 

conjugate and transpose operation of SD multiplied by <3?d- For the second constraint, it means minimizing 

correlation of off-diagonal elements of matrix which is conjugate transpose of SD multiplied by Done 

with it like this, the sensing dictionary W is insensitive to scattering mechanisms. In other words, the 

model mismatch can be imitated greatly. In the following, it illustrates problem P 2 being a convex 
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problem. Noting that W and §d are divided column by column and W = [wi, W2, • • • ,wjy], &d 
[<Pdii 4>d2, ■ ■ ■ j <PdN]- Denoting column vector (f>dh <t>dl an d as follows, respectively, 

,T 



hi = [^dif ^dif] (8) 



(9) 



w z = [K(wO T 9(w,) TlT 



let, 



/i(W) = \\l-diag(W H ^ d )\\ oo -b 1 
Noting that the first inequality constrains in problem P2, 

\\l-diag(W H <f d )\\ oo <b 1 
which is equivalent as the following constraints, 



(10) 



(11) 



(12) 



i-wf 4> d i 

l-W? <h d 2 



J N <PdN oo 



< b 



(13) 



Without loss of generality, it is supposed that the absolute of Zth component is the largest in (12), so 



it has, 



l-wf 

l-wi 



1— w 



NVdN oo 



|1 — w 



Using equation (9), (10) and (11), we can get 



Wl<J>dl + 3 ■ Wj0. 



V// 



Udi\ 



w« [4>di <j) d i] 



(14) 



(15) 



Obviously, it is a second order cone about W. Hence, /i(W) is a convex function about W. For the 
second constraint condition, let 



Obviously, 



f 2 (W) = \\vec ((W^/m) IL " 6 2 (k ± I) 



(16) 



(17) 



Done with similar derivation procedure of /i(W), it is easy to show that /2(W) is also a convex 
function about W and objective function is affine function w.r.t. b\, 62 for a given 7. Hence P2 is convex. 
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4.2. The Proposed Algorithm 

The proposed algorithm (OMP-SD) is also a greedy algorithm but different from OMP and A-OMP 
For OMP-SD, at each iteration, it requires a two-step search to select an atom. First, it determines the 
offset index of atom in SD W, which is not sensitive to scatter mechanism. And then, it further to 
determine the specific scatter mechanism in dictionary <I>. After the two-step procedures, an actual atom 
is picked out. The OMP-SD is described as follows. 

Algorithm 2 :OMP-SD 
Input: 

The measurement vector, y; 
The dictionary, 3>i, $ 2 ■ ■ ■ , $D, W 
the err threshold, e; 
Main Procedures: 

1: To initialize the residual ro = y and initialize the subscribe set cq is empty, set % = 1. 
2: To find the matrix I\ that solves the maximization problem 

i A I H I 

ti = max \ w t r.;_i| 

where, 

r = [$i(to $ 2 (U) ••• $ D (ti)] 

3: To solve the maximization problem 

£ = max Ir^rj-i I , d G A 

d 

and update q = q_i U {£j}. Where £« = (£ — l)iV + tj. 
4: Let Pj = $(cj)($(cj) H< I>(ci))~ 1< I > (cj) i/ denote the projection onto the linear space spanned by the 

elements of <3?(cj). Update = (I — Pj)y. 
5: If the stopping condition is achieved (i.e., ||rj||2 < e), go to 6. Otherwise, set i = i + 1 and go back 

to 2. 

6: Pick out the range scattering cells w.r.t set Cj. 

7: Calculate the scattering intensity in these range cells determined in the previous step with P^y. 
8: Reconstruct SRP using the scattering intensity and range scattering cells. 
9: Return SRP. 
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TABLE I: Geometry parameters for example scattering geometries ifTTl 



Value of ad 


Scatter mechanisms 


-1 


corner diffraction 


-0.5 


edge diffraction 





point diffraction; straight edge specular 


0.5 


singly curved surface reflection 


1 


late plate at broadside; dihedral 



As far as computational complexity is concerned, it requires DN times correlation operators to select 
an atom in OMP whileas it just requires N + D times for the proposed algorithm. So it requires about 
0(M(N + D)K) flops. It is approximate to the simplified model in which requires N times. Similarly 
to in section 3 discussed, once the subscribe set is determined with the proposed algorithm, the SRP can 
be recovered with LS solution, too. 

5. Simulation and Experimental Results 

In this section, 10000 trails Monte Carlo simulation has been done to illustrate the previous discussions. 
Assume the SFR operates at the following condition. Five scattering mechanisms are considered (i.e., 
a.d 6 ^ — { — 1,-0.5,0,0.5,1}). An example scattering geometries and the corresponding scattering 
parameters are shown in Tab. [n] 

In (4), it is assumes that target is stationary in one CPI, and the distance L in Fig. 1 is regarded as 
constant. In simulation, set the tq = and hence, the rath row and nth column element in (4) is rewritten 
as 

[%] m ,n = G p [j(l + mA/// )p • exp(-j2vr/ m n/(AfA/))) (18) 

The range of the measured frequency band is from L band to S band ( i.e., from 1GHz to 4GH ), 
where the start frequency is /o = 1GHz and frequency step size Af = 10MHz. The number of pulses 
M = 300. And it assumes that the target is 5m length. Five scatterers are located on 0.3m, 0.85m, 
2.0m, 3.25m and 4m to target front-end, respectively. All scatterers have same intensity. What's more, it 
assumes that the stationary scatterer centers are present on the grid points. In each measurement, only 30 
returned pulses are measured in one CPI (i.e., 300 pulses). The measurement vector y is contaminated by 
AWGN with SNR = 20dB, 15dB, lOdB, 5dB and noiseless situation, respectively. In order to explain 
the essence of model mismatch, mutual incoherence property (MIP) is introduced which is defined as 
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Fig. 2: IAI in original dictionary 



TABLE II: Time consuming simulation results 



SRP Algorithm 


time(s) 


OMP 


1044 


A-OMP 


109 


OMP-SD 


121 



the same as in article [23, 

^ = II jf ill II (19) 

\\(j)f \\ 2 ■ WfaWz 

Noting that each atom in dictionary is normalized, hence Eq. (19) can be rewritten as another form, 

/i($)= max \4>f(j)j\ (20) 

l^.i,j^.n 

Simulation results are shown in Fig. 2 ~ Fig. 8. and it presents several remarks in the following. 

Remark 1: The IAI between the original dictionary and SD are shown in Fig. 2 and Fig. 3, respectively. 
In the original dictionary, the IAI minimum is 0.0762 and the maximum of inner-atom cross-correlation 
is 0.3872. However, both of them are 0.2248 (i.e., b\ = hi = 0.2248 in (7a)) for the sensing dictionary. 
It improves the autocorrelation between atoms and mitigates model mismatch at the same time. 

Remark 2: Fig. 4 shows that the success recovery probability is a monotonic decreasing relative 
to SNR for the three algorithms (i.e., OMP, A-OMP and OMP-SD). It is easy to understand that 
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Fig. 3: IAI in sensing dictionary 
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Fig. 4: Success recover probability w.r.t SNR 



the OMP has the best recovery performance because it is match model and the A-OMP has worst 
recovery performance because of its model mismatch. However, the proposed method (OMP-SD) has an 
approximate performance compared to OMP and approximate computational complexity to A-OMP and 
it is confirmed in Tab. [n] It should be noted that there is an exception for small SNR (< lOdB). When 
SNR is 5dB in Fig. 4, all of the three algorithms have a lower success recovery probability (less than 30 
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Fig. 5: Cumulative distribute error in noiseless 
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Fig. 6: Cumulative distribute error SNR = 20dB 



percent). Hence, objectively speaking, it is a drawback for these algorithm. But in moderately high SNR 
settings (i.e., greater than 15dB), the proposed algorithm has outstanding performance. 

Remark 3: For the noiseless and three different SNR settings, Fig. 5 ~ Fig. 8 show the cumulative 
distribute errors (i.e., CDE). It is widely used to evaluate recover performance in CS community such 
as ll25l . From Fig. 5 ~ Fig. 8 we can see the match model is best, while the mismatch model is worst 
although it requires least computation amount. However, OMP-SD shows that it has an approximate 
values of CDE compared to OMR However, it has to point out that all of the three algorithms are not 
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Fig. 7: Cumulative distribute error SNR = 15dB 
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Fig. 8: Cumulative distribute error SNR = lOdB 



suitable for low SNR (< lOdB) settings. 

Remark 4: In Tab. |n| for the same computer platform, there are sums of 10000 Monte Carlo trails 
time consuming results for OMP, A-OMP and OMP-SD, respectively. It confirms that OMP has the most 
computational cost but it is approximately computational cost between A-OMP and OMP-SD. Both of 
them have much lower computation cost compared with OMP. 

6. Conclusion and Future Work 

In this paper, a fast algorithm to synthesize range profile is proposed. For the SFR system in GTD 
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model, the HRRP synthesis can be converted to solve a sparse approximation problem over redundant 
dictionaries. Different from A-OMP, the model mismatch is mitigated with SD. Better than OMP, the 
computational complexity is reduced. Finally, simulation results show the proposed algorithm is valid for 
both noiseless and noisy settings. 

However, it just presents the fast algorithm to recover parameters in GTD model. In the future work, 
it will derive the ERC. Secondly, in this paper, it just presents how to estimate the support of sparse 
vector. For the scatterer intensity in GTD model, it simply exploits a plain LS method to recover it. In 
fact, there are a great many algorithms to estimate its magnitude such as biased estimation techniques 
E71 . etc. All of them will be considered in our future work. 
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